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ABSTRACT 

Context. Particle trapping in local or global pressure maxima in protoplanetary disks is one of the new paradigms in the theory of 
the first stages of planet formation. However, finding observational evidence for this eff'ect is not easy. Recent work suggests that the 
large ring-shaped outer disks observed in transition disk sources may in fact be lopsided and constitute large banana-shaped vortices. 
Aims. We wish to investigate how eff'ective dust can accumulate along the azimuthal direction. We also want to find out if the size- 
sorting resulting from this can produce a detectable signatures at millimeter wavelengths. 

Methods. To keep the numerical cost under control we develop a 1 + lD method in which the azimuthal variations are treated sepa- 
rately from the radial ones. The azimuthal structure is calculated analytically for a steady- state between mixing and azimuthal drift. 
We derive equilibration time scales and compare the analytical solutions to time-dependent numerical simulations. 
Results. We find that weak, but long-lived azimuthal density gradients in the gas can induce very strong azimuthal accumulations of 
dust. The strength of the accumulations depends on the Peclet number, which is the relative importance of advection and diff'usion. 
We apply our model to transition disks and our simulated observations show that this eff'ect would be easily observable with ALMA 
and in principle allows to put constraints on the strength of turbulence and the local gas density. 

Key words, accretion, accretion disks - protoplanetary disks - stars: pre-main- sequence, circumstellar matter - planets and satellites: 
formation - submillimeter: planetary systems 



1. Introduction 

The process of planet formation is thought to start with the 
growth of dust aggregates in a protoplanetary disk through coag- 
ulation (see e.g., the review by Blum & Wurm 2008). The idea is 
that these aggregates get successively bigger until either gravo- 
turbulent processes set in (Goldreich & Ward 1973; Johansen 
et al. 2007; Cuzzi et al. 2008) or planetesimals are formed di- 
rectly via coagulation (Weidenschilling 1977; Okuzumi et al. 
2012; Windmark et al. 2012). One of the main unsolved prob- 
lems in these scenarios is the problem of excessive radial drift 
(Weidenschilling 1977; Nakagawa et al. 1986; Brauer et al. 
2007). The origin of this problem lies in the sub-Keplerian mo- 
tion of the gas in the disk, caused by the inward pointing pres- 
sure gradient. The dust particles in the disk, however, feel only 
the friction with the sub-Keplerian gas and as they grow to 
larger sizes, the reduced surface-to-mas s ratio causes them to 
fall toward the star with speeds up to 50 m s"^ (Whipple 1972; 
Weidenschilling 1977). Thus, if a particle grows, it will sooner or 
later get "flushed" toward the star before it can grow very large 
(Brauer et al. 2008; Birnstiel et al. 2010, 2012). This is what we 
call the radial drift barrier, and this still poses one of the main 
unsolved problems of the early phases of planet formation. 

A possible solution to this problem might lie in the con- 
cept of "particle traps". This idea has been proposed already 
some time ago in the context of anticy clonic vortices by Barge 
& Sommeria (1995) and Klahr & Henning (1997), as well as 
in/around turbulent eddies (Johansen & Klahr 2005; Cuzzi & 
Hogan 2003). These vortex or eddy-related local pressure max- 
ima act as particle traps because particles tend to drift always to- 



wards higher pressure. In the case of vortices and eddies they act 
on small scales. Particle traps on a global disk scale have been 
proposed as well (Whipple 1972; Rice et al. 2006; Alexander & 
Armitage 2007; Garaud 2007; Kretke & Lin 2007; Dzyurkevich 
et al. 2010) and shown to be conducive to planet formation 
(Brauer, Henning, & Dullemond 2008). Intermediate scale par- 
ticle traps may also occur from magnetorotationally driven tur- 
bulence: the so-called zonal flows (Johansen et al. 2009). 

If we want to observationally test whether this scenario of 
particle trapping actually occurs in nature, we are faced with a 
problem. The Earth-forming region around a pre-main sequence 
star is usually too small on the sky to be spatially resolved suffi- 
ciently well to test this trapping scenario. Moreover, the optical 
depth of this inner disk region is likely to be too large to be 
able to probe the mid-plane region of the disk. Fortunately, what 
constitutes the "meter size drift barrier" at 1 AU is a "cm size 
drift barrier" at ~ 50 AU. Those disk regions are optically thin at 
millimeter (mm) wavelength and particles in the mm size range 
can be spectroscopically identified by studying the mm spectral 
slope (Testi et al. 2001; Natta et al. 2004; Ricci et al. 2010a,b). 
So the goal that has been pursued recently is to identify obser- 
vational signatures of dust particle trapping of millimeter- sized 
particles in the outer regions of disks, as a proxy of what hap- 
pens in the unobservable inner regions of the disk (Pinilla et al. 
2012b,a). The Pinilla et al. (2012a) paper suggests that the huge 
mm continuum rings observed in most of the transition disks 
(Pietu et al. 2006, Brown et al. 2008, Hughes et al. 2009, Isefla 
et al. 2010, Andrews et al. 2011), may in fact be large global 
particle traps caused by the pressure bump resulting from, for 
example, a massive planet opening up a gap. 
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We have focused in these papers on the intermediate scale 
(zonal-flow-type) and the the large scale (global) pressure 
bumps, simply because current capabilities of mm observatories 
(including ALMA) are not yet able to resolve small scale struc- 
tures such as vortices. The vortex trapping scenario thus appears 
to remain observationally out of reach. However, a closer look 
at the mm maps of transition disks suggest that some of them 
may exhibit a deviation from axial symmetry. For instance, the 
observations presented in Mayama et al. (2012) or the mm im- 
ages of Brown et al. (2009) suggests a banana shaped lopsided 
ring instead of a circular ring. Regaly et al. (2012) proposed that 
these banana- shaped rings are in fact a natural consequence of 
mass piled up at some obstacle in the disk. Once the resulting 
ring becomes massive enough, it becomes Rossby-unstable and 
a large banana- shaped vortex is formed that periodically fades 
and re-forms with a maximum azimuthal gas density contrast 
of a few. Regaly et al. (2012) showed that this naturally leads 
to lopsided rings seen in mm wavelength maps (see also ear- 
lier work by Wolf & Klahr 2002) on radiative transfer predic- 
tions of observability of vortices with ALMA). The formation of 
such Rossby-wave induced vortices was demonstrated before by 
Lyra et al. (2009), who showed that they may in fact (when they 
are situated much further inward, in the planet forming region) 
lead to the rapid formation of planetary embryos of Mars mass. 
Sandor et al. (2011) subsequently showed that this scenario may 
rapidly produce a 10 Earth mass planetary core. 

The goal of the current letter is to combine the scenario of 
forming a lopsided gas ring (e.g., Regaly et al. 2012) with the 
scenario of particle trapping and growth presented by Pinilla 
et al. (2012a). In Section 2, we will outline the physical ef- 
fects involved and derive analytical solutions to the dust distribu- 
tion along the non-axisymmetric pressure bump and Section 3, 
we will test the observability of these structures in resolved 
(sub-)mm imaging and in the mm spectral index. Our findings 
will be summarized in Sect. 4. 



2. Analytical model 

Dust particles embedded in a gaseous disk feel drag forces if they 
move relatively to the gas. The radial and azimuthal equations 
of motion have been solved for example by Weidenschilling 
(1977) or Nakagawa et al. (1986) for the case of a axisymmet- 
ric, laminar disk. It was found that particles drift inward to- 
wards higher pressure. In this paper, we will focus on the case 
where a non-axisymmetric structure has formed a long lived, 
non-axisymmetric pressure maximum in the disk. This pressure 
maximum is able to trap inward- spiralling dust particles. Like in 
the aforementioned works we can solve for a stationary drift ve- 
locity, but unlike in Nakagawa et al. (1986), the radial pressure 
gradient is zero at the pressure maximum while the azimuthal 
pressure gradient can be diff'erent from zero. This leads again to a 
systematic drift motion of the dust particles towards the pressure 
maximum, but now in azimuthal instead of in radial direction. 

The equations of motion in polar coordinates relative to the 
Keplerian motion become (see Nakagawa et al. 1986) 



-Apg(Ud,r - Ug^r) + 2Qk%0 (1) 
-Apg(Wd,0 - Wg,^) - ^^kUd,r (2) 
= -Apd(%r - Ud,r) + 2Qk%0 (3) 



dt 



- -Apd(wg,0 - Wd,(^ 



1 dP 



(4) 



where t/d,r and t/d,0 are the r and components of the dust veloc- 
ity, respectively, u^^r and Wg,^ the r and components of the gas, 
Qk the Keperian frequency, P the gas pressure and pd and pg 
the dust and gas densities. A denotes the drag coefficient (see, 
Nakagawa et al. 1986, equations 2.3 and 2.4). Solving above 
equations for the velocity along the direction at the mid-plane 
{z = 0) gives 



1 dP 



St -h St"^ (1 -h X)2 PgVk#' 



(5) 



where X = Pd/Pg is the dust-to-gas ratio, Vk the Keplerian ve- 
locity, and the Stokes number^ is given by 



St: 



Psa 



:nk. 



(6) 



with ps as internal density of the dust, particle radius a, and the 
isothermal sound speed Cs. 

Thus, dust is advected with the velocity given in Eq. 5 but 
it is also turbulently stirred. Together, the evolution of the dust 
density along the ring is then described by 

^ = |;(pd«d,,)-|;(i)pg|;(^)), (7) 

where y = r0 is the coordinate along the ring circumference. 
We use a dust diffusion coefficient D according to Youdin & 
Lithwick (2007), D = Dg^J(l -h St^), where we assume the 
gas diffusivity to be equals to the gas viscosity, taken to be 
y = QftCg/Qk, with at as the turbulence parameter (see Shakura 
& Sunyaev 1973). Equation 7 can be integrated forward in time 
numerically, but assuming that the turbulent mixing and the drift 
term have reached an equilibrium and also assuming a low dust- 
to-gas ratio, we can analytically solve for the dust density in a 
steady-state between mixing and drifting, which yields 



Pd(y) = C-pg(3;)-exp 



St(3;) 



at 



(8) 



where C is a normalization constant and St(y) is the Stokes num- 
ber which depends on y via the changes in gas density. Eq. 8 thus 
predicts the distribution of dust for any given profile of the gas 
density pg. The contrast between the position of the azimuthal 
pressure maximum and its surrounding then gives 



^min 



exp 



3^rmn _ 3^0 



at 



(9) 



which is plotted in Fig. 1. St"''''' and St'"''' are the Stokes num- 
bers at the pressure maximum and minimum, respectively (note: 
3^max ^ stmin^j ^^^^ ^^^^ ^^^^ particlcs Stokcs num- 

ber becomes larger than the turbulence parameter at, the dust 
concentration becomes much stronger than the gas concentra- 
tion. However, the time scales on which these concentrations are 
reached can be significant. In order to get an estimate for this 
time scale, we compare the advection time scale ^adv = L/u and 
the diffusion time scale ^diff = L^/D, with velocity u and length 



^ St =1 typically corresponds to particles of mm to cm sizes in the 
outer disk. For typical disk conditions a ^ 0.4 cm • St • Eg/1 g cm~^. 
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Fig. 1. Contrast between the dust density in the azimuthal max- 
imum and its surroundings for different gas density contrasts 
^max^pmm^^ Stokes numbers St, as derived in Eq. 9. Turbulence 

parameters are = 10"^ (soHd), and = 10"^ (dashed). 



Fig. 2. Azimuthal steady state solutions for small (red line) and 
large dust grains (black line) for a given sinusoidal gas profile 
(black line, identical with red line). Dashed lines show numerical 
solutions at 1,2, and 5 advection time scales. 



scale L. The ratio of these time scales, known as the Peclet num- 
ber Pe, which, using Eq. 5, can be written as 



Pe 



^diff ^ St L ^Pg 
^adv o:t Pg dy 



StApg 

at Pg 



(10) 



It describes the relative importance of advection and diff'usion 
and confirms that dust accumulations occur only for particles 
with St > at, because otherwise, diff'usion dominates over ad- 
vection which means that variations in the dust-to-gas ratio are 
being smeared out. It also shows that for those large particles 
the advection time scale is the shorter one, thus setting the time 
scale of the concentration process, which can be written as 



^adv 



TT^ ^HJ^ 1 



(11) 



^St \ r / Qk 

where H = c^/fl^ is the pressure scale height and we have used 
a mean velocity u = ud(f>\ = ^^St, with 



^St = |arccot(St"^^^) - arccot (St"^"^)| 



(12) 



which for St < 1 simplifies to ^St = St^^ - St^^^ Further, we 
define the pressure maximum and minimum to be at = and 
(p = 71, respectively. As an example, at 35 AU, for H/r = 0.07, a 
Stokes number of 0.2 and a gas density contrast of p^^^/p^^^ = 

2, the time scale is 3 x 10^ years, but could be as short as 10^ 
orbits for optimal conditions. Any gas structure therefore has to 
be long-lived to cause strong asymmetries in the dust, making 
e.g. asymmetries caused by a planet or long lived vortices the 
best candidates (Meheut et al. 2012). If such an accumulation is 
formed and the gas asymmetry disappears, it still takes ^diff to 
"remove" it, which at 35 AU is of the order of Myrs. It remains 
to be shown whether short-lived, but reoccurring structures like 
zonal flows are able to induce strong dust accumulations. 

3. Simulated Observations 

In the following, we will evaluate whether the dust structures we 
expect would be observable in resolved mm images of ALMA. 



As radial gas surface density profile ^g(r) and temperature T(r) 
we use the results presented in Pinilla et al. (2012a). These sim- 
ulations represent a disk of mass M^isk = 0-05 M© around a solar 
mass star with a 15 Jupiter-mass planet at 20 AU. For simplicity 
the two-dimensional gas surface density is then taken to be 



i:g(r,0) = i:g(r). 



1 -h A(r) ' sin 



in(0-f) 



c-h 1 



• exp 



(r-Rsf 



2 if 2 



(13) 
(14) 



where c = I^g,max/^g,min is the largest contrast of the gas surface 
density, taken to be 1.5. is the position of the radial pressure 
bump. The dust size distribution Ed(r, a) was also taken from the 
simulations of Pinilla et al. (2012a), and distributed azimuthally 
using the analytical solution from Eq. 8 (see Fig. 2). We also con- 
firmed the analytical solution and time scales by solving Eq. 7 
numerically, as shown in Fig. 2. This analytical solution strictly 
only holds at the position of the radial pressure bump, but since 
most of the mm emission comes from the large grains which in 
the simulations of Pinilla et al. (2012a) are trapped near the ra- 
dial pressure maximum, this should be a reasonable approxima- 
tion. Full two-dimensional simulations will be needed to confirm 
this and to investigate the eff'ects of shear. 

To compare directly with current ALMA observations, we 
calculate the opacities for each grain size at diff'erent wave- 
lengths and assume spherical silicate grains with optical con- 
stants for magnesium-iron grains from the Jena database^. The 
continuum intensity maps are calculated assuming that in the 
sub-mm regime, most of the disk mass is concentrated in the 
optically thin region. We assume the same stellar parameters as 
in Pinilla et al. (2012a), azimuthally constant temperature T(r), 
typical source distances (d = 140 pc) and zero disk inclination. 
We run ALMA simulations using CAS A (v. 3.4.0) at 345 GHz 
(band 7) and 675 GHz (band 9), shown in Fig. 3. We consider 
2 hours of observation, the most extended configuration that is 
currently available with Cycle 1 , generic values for thermal and 
atmospheric noises and a bandwidth of Ay = 7.5 GHz for con- 
tinuum. At these two diff'erent frequencies, it is possible to detect 



http ://w w w. astro .uni-j ena. de/Laboratory/Database/databases .html 
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Relative J2000 Right Ascension (arcsec) 

Fig. 3. ALMA simulated images at 345 GHz with an observation 
time of 2 hours. The total flux of the source is 0.13 Jy and the 
contour lines are at {2, 4, 6, 8} the rms (cr = 0.22 mJy). 

and resolve regions where the dust is trapped creating a strong 
azimuthal intensity variation. 

The spectral slope CKmm of the spectral energy distribution 
Fy oc is directly related to the dust opacity index at these 
long wavelengths (e.g. Testi et al. 2003), and it is interpreted in 
terms of the grain size (amm ^ 3 implies mm sized grains). 
With the simulated images of Fig. 3, we compute the aj^m map 
(Fig. 4), considering an antenna configuration that provides a 
similar resolution of ~ 0.16" (~ 22AU at 140 pc) for each band. 
This resolution is enough to detect aj^m variations along the az- 
imuth, confirming that those are regions where dust accumulates 
and grows due to the presence of an azimuthal pressure bump. 

4. Summary and conclusions 

We have shown that weak, but long-lived azimuthal asymmetries 
in the gas density can cause very efficient accumulation of dust 
at the position of the azimuthal pressure maximum. We have de- 
rived analytical steady-state solutions for the dust distribution for 
any given azimuthal gas density distribution and the time scales 
on which these distributions develop. Good agreement has been 
found between the solutions and numerical simulations. 

For this dust concentration mechanism to work, particles 
must have grown to larger sizes (St > at) such that the azimuthal 
drift becomes stronger than the turbulent diff'usion. This typi- 
cally corresponds to particles of sub-mm to cm sizes. The strong 
concentration of the largest grains leads to a size-sorting which 
is observable via low spectral indices at mm wavelengths and 
also in lopsided banana shaped structures in resolved (sub-)mm 
images. Finding the size range where the bifurcation between 
concentration and diff'usion happens would put constraints on 
the turbulence strength and the local gas density of the disk. 

Acknowledgements. We thank Ewine van Dishoeck, Simon Bruderer, Nienke 
van der Marel, Geoffrey Mathews, Hui Li and the referee for useful discussions. 



References 

Alexander, R. D. & Armitage, P. J. 2007, MNRAS, 375, 500 
Andrews, S. M., Wilner, D. J., Espaillat, C, et al. 2011, ApJ, 732, 42 
Barge, R & Sommeria, J. 1995, A&A, 295, LI 




_^ , , : 1 , , , , 1 2 

0.8 0.6 0.4 0.2 -0.2 -0.4 -0.6 -0.8 



Relative J2000 Right Ascension (arcsec) 

Fig. 4. Spectral index map Qfmm using simulated images at band 7 
and 9. The antenna configuration is chosen such that the angular 
resolution is similar for both bands ~ 0.16" (~ 22AU at 140 pc). 

Birnstiel, T., Dullemond, C. R, & Brauer, R 2010, A&A, 513, 79 
Birnstiel, T., Klahr, H., & Ercolano, B. 2012, A&A, 539, 148 
Blum, J. & Wurm, G. 2008, ARA&A, 46, 21 

Brauer, R, Dullemond, C. R, Johansen, A., et al. 2007, A&A, 469, 1169 
Brauer, R, Henning, T., & Dullemond, C. R 2008, A&A, 487, LI 
Brown, J. M., Blake, G. A., Qi, C, Dullemond, C. R, & Wilner, D. J. 2008, ApJ, 
675, L109 

Brown, J. M., Blake, G. A., Qi, C, et al. 2009, ApJ, 704, 496 
Cuzzi, J. N. & Hogan, R. C. 2003, Icarus, 164, 127 
Cuzzi, J. N., Hogan, R. C, & Shariff, K. 2008, ApJ, 687, 1432 
Dzyurkevich, N., Flock, M., Turner, N. J., Klahr, H., & Henning, T. 2010, A&A, 
515,70 

Garaud, R 2007, ApJ, 671, 2091 
Goldreich, R & Ward, W. R. 1973, ApJ, 183, 1051 
Hughes, A. M., Andrews, S. M., Espaillat, C, et al. 2009, ApJ, 698, 131 
Isella, A., Natta, A., Wilner, D., Carpenter, J. M., & Testi, L. 2010, ApJ, 725, 
1735 

Johansen, A. & Klahr, H. 2005, ApJ, 634, 1353 

Johansen, A., Oishi, J. S., Low, M.-M. M., et al. 2007, Nature, 448, 1022 

Johansen, A., Youdin, A., & Klahr, H. 2009, ApJ, 697, 1269 

Klahr, H. H. & Henning, T. 1997, Icarus, 128, 213 

Kretke, K. A. & Lin, D. N. C. 2007, ApJ, 664, L55 

Lyra, W., Johansen, A., Klahr, H., & Piskunov, N. 2009, A&A, 493, 1125 

Mayama, S., Hashimoto, J., Muto, T., et al. 2012, ApJL, 760, L26 

Meheut, H., Keppens, R., Casse, R, & Benz, W. 2012, Astronomy & 

Astrophysics, 542, 9 
Nakagawa, Y., Sekiya, M., & Hayashi, C. 1986, Icarus, 67, 375 
Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 2004, A&A, 416, 

179 

Okuzumi, S., Tanaka, H., Kobayashi, H., & Wada, K. 2012, ApJ, 752, 106 
Pietu, v., Dutrey, A., Guilloteau, S., Chapillon, E., & Pety, J. 2006, A&A, 460, 
L43 

Pinilla, R, Benisty, M., & Birnstiel, T. 2012a, A&A, 545, 81 
Pinilla, R, Birnstiel, T., Ricci, L., et al. 2012b, A&A, 538, 114 
Regaly, Z., Juhasz, A., Sandor, Z., & Dullemond, C. R 2012, MNRAS, 419, 1701 
Ricci, L., Testi, L., Natta, A., & Brooks, K. J. 2010a, A&A, 521, 66 
Ricci, L., Testi, L., Natta, A., et al. 2010b, A&A, 512, 15 
Rice, W. K. M., Armitage, R J., Wood, K., & Lodato, G. 2006, MNRAS, 373, 
1619 

Sandor, Z., Lyra, W, & Dullemond, C. R 2011, ApJL, 728, L9 

Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 

Testi, L., Natta, A., Shepherd, D. S., & Wilner, D. J. 2001, ApJ, 554, 1087 

Testi, L., Natta, A., Shepherd, D. S., & Wilner, D. J. 2003, A&A, 403, 323 

Weidenschilling, S. J. 1977, MNRAS, 180, 57 

Whipple, F. L. 1972, From Plasma to Planet, 211 

Windmark, F., Birnstiel, T., Ormel, C. W, & Dullemond, C. R 2012, A&A, 544, 
L16 

Wolf, S. & Klahr, H. 2002, ApJ, 578, L79 

Youdin, A. N. & Lithwick, Y. 2007, Icarus, 192, 588 



4 



